Multi-plectoneme phase of double-stranded DNA under torsion 
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Abstract 

We use the worm-like chain model to study supercoiling of DNA under tension and torque. The 
model reproduces experimental data for a much broader range of forces, salt concentrations and 
contour lengths than previous approaches. Our theory shows, for the first time, how the behavior 
of the system is controlled by a multi-plectoneme phase in a wide range of parameters. This 
phase does not only affect turn-extension curves but also leads to a non-constant torque in the 
plectonemic phase. Shortcomings from previous models and inconsistencies between experimental 
data are resolved in our theory without the need of adjustable parameters. 



The DNA contained in every cell of all higher organisms is hundred times longer than 
the cell diameter: to fit inside it has to fold. This is a challenging problem since DNA is 
a semi-flexible polymer, making it hard to confine it in small spaces. On the other hand, 
local unfolding of DNA has to be efficient, as it plays a key role in the transcription and 
replication of the genome. Unfolding is achieved by stretching and twisting the molecule: 
unraveling how DNA reacts to them is crucial to understand cellular activities. 

The relevant mechanical properties of DNA have been studied with single molecule tech- 
niques, where individual molecules are stretched and/or twisted under physiological condi- 
tions. The stretching and bending elasticity, in the absence of twisting, has been investigated 
through measurements of the force-extension relation of DNA [ ] and theories based on the 
worm-like chain (WLC) model have successfully explained the experimental results [2] [3]. 
The WLC model [i] is a coarse-grained approximation where the particular sequence of 
basepairs (bp) is hidden by treating the DNA as a homogeneous semifiexible polymer. 



To understand stretched DNA under torsional stress [■^>-7] several models based on the 
WLC framework were proposed. However, they were either purely mechanical [ ], involved 
non-linear elasticity [9], phenomenological [I!)], aimed only at a certain region of the experi- 
mental data [ ] or had to invoke a non-canonically reduction of the DNA charge density [12] . 
The outcome of the experiments still remains poorly understood. 

In this Letter we present a theory, based on the WLC model, without any of the afore- 
mentioned shortcomings. The results describe experimental data accurately (see Fig. 2 for 
an example). Up to now it has been assumed that under high twist DNA reduces its torque 
through the formation of a single plectoneme, see Fig. 1. We show here for the first time 
how thermal fluctuations lead to a multiple plectoneme phase instead. We demonstrate its 




Figure 1. Geometry of the plectoneme. 
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impact on the turn-extension slope, see Fig. 3, and on the torque, see Fig. 5. 
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Figure 2. The turn-extension plots of a 600 nm DNA chain at 320 mM ionic strength for various 
tensions between 0.25 and 4pN. Comparison between theory and experiments. Experimental data 
from [12]. 
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Figure 3. The results of the theory with and without the possibility to form more than one 
plectoneme are presented alongside the experimental results (3pN, 20 mM, Experimental data 
from [12]. 



In the experiments the DNA is anchored at one end to a surface and at the other end 
to a magnetic [ ] or optical [ ] tweezer. This allows to control the tension and the linking 
number (n), the number of turns inserted, at the same time. Increasing n at constant tension 
yields turn-extension plots like the ones shown in Fig. 2. Initially most of n goes into twist 
{Tw) of the molecule while the end-to-end distance remains approximately constant. Then 
a transition — dependent on the salt concentration, the DNA length and the applied force 
and often accompanied by a jump [u] — signals the formation of a plectoneme. From this 



point onward all the additional n is stored inside the growing plectoneme as writhe {Wr), a 
quantity related to the path of the polymer. Writhe, twist and linking number are connected 
through White's relation [13], n = Tw + Wr. 

We assume that the legs and the end loop form a homoclinic solution, characterized by 
the parameter t G [0,1], as described in [1 []. t = corresponds to a straight rod, and 
t = 1 to the homoclinic loop. Inside the homoclinic solutions, at the point of (non-zero) 
minimum distance of symmetric points, we insert a plectoneme. Such a minimum exists in 
the range 0.804 > t > 1; half this distance sets the plectoneme radius R{t). The homoclinic 
solution stores some fixed amount of writhe, Wri = 2(arcsint)/7r. Moreover its bending 
and potential energy sum up to Ei = 8FXt. Here F is the tension, Pt = A/ksT the DNA 
bending persistence length and A = \/ AjF . 

On the other hand, the plectoneme has a bending energy density e?, = Aco-^^ ajlF^it) 
and a contribution to the potential energy density of / where a denotes the angle of the plec- 
toneme. Fig. 1. When properly accounting for the presence of the end loop, the plectoneme 
writhe density is exactly given by wr^if) = sin a cos a/ 271 R{t). The electrostatic interaction 
in the plectoneme has a free energy density /ei,o described by Ubbink and Odijk [15]. The 
effective charge density in this contribution is calculated on the basis of a charge density of 
2 charges/0.34 nm with the method described in [16]. 

The zero-temperature energy density of a DNA chain with contour Lc containing m 
plectonemes of total relative length Ip = Lp/L^ is given by 



where the last term is the twist contribution to the energy, and Pc is the torsional persistence 
length. However, to properly account for the experimental situation the zero-temperature 
analysis is not sufficient. 

Thermal fluctuations lead to three different contributions: the first affects the DNA legs, 
the second acts within plectonemes and the last is related to the number of plectonemes and 
their position and length distribution. Outside the plectonemes, thermal fiuctuations modify 
the shape of the DNA; for a given torque r they shorten the DNA end-to-end distance by 
a factor p(F, r) [ ] and give origin to thermal writhe that lowers the twist energy density, 
which can be expressed as a renormalized torsional persistence length P^[X) [17]. 



eo(m, Ip) = IpCb + m-^ + lp{F + /ei,o) 




(1) 
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Inside the plectonemes, following Ref. [1^)] we consider fluctuations in two directions. We 
denote by a,,, {ap) the standard deviation in the radial (pitch) direction. The fluctuations 
in the pitch direction are set by the geometry: dp = Tripsin a, see Fig. 1, while in the 
other direction they are set by the electrostatic repulsion. Their contribution changes the 
electrostatic energy, up to flrst order, to fei = fe\,o x exp(4/t^(T^) [15], where is the Debije 
length. The presence of twist couples the two directions of the fluctuations, an effect not 
studied before [18]. This result in a new effective deflection length: 

(A. + Ap)(A2 + A2) 

where Ar.,p = (P6(o"^p))^''^ are the deflection lengths of confinement as defined in Ref. [19]. 
The resulting twist energy density, renormalized with P^^{X), depends on the confinement. 
Since twist diffusion happens on a very short time scale, the twist energy density in the legs 
and in the plectonemes should be the same. This non-trivially couples the linking number 
density between legs and plectonemes. 

As a consequence of thermal fluctuations, the plectonemes' path is shortened by Pstr? its 
bending energy density (cf,) and bare writhe density (wrb) are modifled; the conflnement of 
the polymer contributes an additional /conf to the free energy [20]. Together with the change 
from /cLo to /el, these modify eo{m,lp) to f{m,lp). 

In the inflnite chain limit, as long as the number of plectonemes stays small, the plec- 
toneme parameters R, a and become independent of Ip since end loop contributions drop 
out. 

The flnal contribution of the fluctuations comes from two combinatorial factors in the 
partition function. They arise from the number of ways the total length of the plectonemes 
can be distributed between them, and the number of ways the plectonemes can be placed 
along the DNA. Since the quantities involved are continuous we need to impose a ^-cutoff 
in our calculations which we assume in the following to be given by the DNA helical repeat, 
^ = 3.4 nm. Assuming hard-core interactions between the plectonemes, this results in the 
partition function 

Z = Zo + Y, I dLpG(m, Lp)e-''^(™''^)^^ (3) 

m=l 



{p{Lc - Lp) - mLioop)'" L"^ ^ 



I ^m. — 1 / — 1 \ I * ^ 
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Here Zq is the partition function when m = and Lioop is the length of a single end loop. The 
first factor of G is the number of ways one can arrange m plectonemes along the DNA; the 
second factor is the number of ways the length Lp can be distributed between m plectonemes. 

When 9„mLioop ^ dnLp the system is in the single-plectoneme state. On the other 
hand, for (9„mLioop ~ dnLp the physics of increasing n cannot be described by the notion 
of plectoneme length growth alone, but a full multi-plectoneme approach is needed. To 
characterize these two states we introduce the multi-plectoneme parameter rj as the difference 
between the writhe efficiencies of loops and plectonemes 

El fp 



V = - (5) 

Wri wTp 

where fp is the free energy density difference between plectoneme and legs, rj is important 
because it enters exponentially the multi-plectoneme parameter C, 

2 

(6) 



V WTp 



One can show that C ^ 1 corresponds to a single-plectoneme state, whereas ^ ^ 1 signals 
the multi-plectoneme phase. In the inset of Fig 4, ^ is displayed as a function of salt 
concentration for different tensions. 
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Figure 4. Phase diagram of the average number of plectonemes as a function of tension and salt 
concentration for a 7.2 jjLvn long chain. Note the shift of the maximum from low tension at high 
salt to high tension at low salt. The inset shows the multi-plectoneme parameter versus salt 
concentration for IpN (blue), 2pN (green) and 3pN (red). 



Theories without the possibility of multiple plectonemes, typically predict, for low salt 
concentrations, a too steep slope of the linear part of the turn-extension curves. As can 
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be seen in Fig. 3 the multi-plectoneme phase accurately describes the experiments, even for 
very low salt concentrations (20 mM, data from [12]). 
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Figure 5. The supercoiUng density-torque plots of a 5600 nm DNA chain at 100 mM ionic strength 
for 0.91 pN, 1.8 pN and 3pN tension. Comparison between theory and torques inferred from ex- 
tension measurements. Data taken from [7] 




Figure 6. The supercoihng density-torque plots of a 725 nm DNA chain at 100 mM ionic strength for 
1 pN, 2 pN and 3 pN tension. Comparison between theory and torques measured using a specially 
crafted cylinder. Data taken from [G] 

Besides the slope, the multi-plectoneme phase influences the torque after the transition. 
In fact, after the transition into the single-plectoneme phase, n is transferred at a fixed 
rate into Ip. This results in a small bump in the torque at the transition, caused by the 
use of a number of turns clamp, followed by a constant plectoneme torque. On the other 
hand, in the multi-plectoneme phase the entropic contribution to the free energy ceases to 
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be linear in n. This explains the difference between torques measured in optical tweezer 
experiments [6] and calculated using Maxwell relations in a magnetic tweezer setup [7] . The 
latter method is based on the assumption of a constant torque after the transition. However, 
for the multi-plectoneme phase our theory predicts a non-constant torque. In Fig. 5 we show 
what our model predicts for the data presented in [/]. To facilitate comparison with the 
original paper, not the linking number, but the supercoiling density is used, which is defined 
as the ratio of the linking number density to the linking density of the two strands of free 
DNA. As can be seen in Fig. 5, the assumption of constant torque underestimates the torque 
difference between the high and low tension curves. The multi-plectoneme phase prediction, 
however, correctly reproduces the torque directly measured in [U], as is shown in Fig. 6. 

A final consequence of the multi-plectoneme phase is the change in the dynamics of 
plectonemes in a chain torsionally loaded. Twist mediated plectoneme length redistribution 
over the plectonemes makes a fast diffusion of plectonemes possible also in the crowded 
environment of the plasmoid in bacteria or a dense chromatin fiber in eukaryotes. The 
implications for cellular processes from transcription to segregation are significant. 

To conclude: the results of the theory show how the multi-plectoneme phase is crucial to 
understand the static and dynamic behavior of DNA under tension and torque. The torque 
in fact is much higher than what was previously computed [~], a fact that could be crucial 
in the understanding of the life processes in which DNA is involved. 

The authors thank Ralf Seidel for providing us with experimental data and Theo Odijk 
and Ralf Seidel for fruitful discussions. We acknowledge discussions with Cees Dekker and 
Marijn van Loenhout. 
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